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Abstract. The techniques employed to solve the interaction of a detector and a 
quantum field typically require perturbation methods. We introduce mathematical 
techniques to solve the time evolution of an arbitrary number of detectors interacting 
with a quantum field. Our techniques apply to harmonic oscillator detectors and can 
be generalized to treat detectors modeled by quantum fields. Since the interaction 
Hamiltonian we introduce is quadratic in creation and annihilation operators, we are 
able to draw from continuous variable techniques commonly employed in quantum 
optics. 



PACS numbers: 03.67.-a, 04.62.+V, 42.50.Xa, 42.50.Dv 
1. Introduction 

The field of quantum information aims at understanding how to store, process, transmit, 
and read information efficiently exploiting quantum resources [1]. In the standard 
quantum information scenarios observers may share entangled states, employ quantum 
channels, quantum operations, classical resources and perhaps more advanced devices 
such as quantum memories and quantum computers to achieve their goals. In order 
to implement any quantum information protocol, all parties must be able to locally 
manipulate the resources and systems which are being employed. Although quantum 
information has been enormously successful at introducing novel and efficient ways of 
processing information, it still remains an open question whether relativistic effects 
can be used to enhance current quantum technologies and give rise to new relativistic 
quantum protocols. 

The novel and exciting field of relativistic quantum information has recently gained 
increasing attention within the scientific community. An important aim of this field 
is to understand how the state of motion of an observer and gravity affects quantum 
information tasks. For a review on developments in this direction see [2]. Recent work 
has focussed on developing mathematical techniques to describe localized quantum 
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fields to be used in future relativistic quantum technologies. The systems under 
investigation include fields confined in moving cavities [3] and wave-packets [4] [5]. 
Moving cavities in spacetime can be used to generate observable amounts of bipartite 
and multipartite entanglement [6,7]. Interestingly, it was shown that the relativistic 
motion of these systems can be used to implement quantum gates [8], thus bridging 
the gap between relativistic-induced effects and quantum information processing. In 
particular, references [7, 8] employed covariance matrix formalism within the framework 
of continuous variables and showed that most of the gates necessary for universal 
quantum computation could be obtained by simply moving the cavity through especially 
tailored trajectories [9]. 

A third local system that has been considered for relativistic quantum information 
processing is the well known Unruh-deWitt detector [10], a point-like quantum system 
which follows a classical trajectory in spacetime and interacts locally with a global 
free quantum field. Such system has been employed with different degrees of success 
in a variety of scenarios, such as in the work unveiling the celebrated Unruh effect 
[10] or to extract entanglement from the vacuum state of a bosonic field. Unruh De- 
Witt detectors seem convenient for relativistic quantum information processing however, 
the mathematical techniques involved, namely perturbation theory, become extremely 
difficult to handel even for simple quantum information tasks such as teleportation 
[11]. The main aim of our research program is to develop detector models which are 
mathematically simpler to treat so they can be used in relativistic quantum information 
tasks. A first step in this direction was taken in [12] where a model to treat analytically a 
finite number of harmonic oscillator detectors interacting with a finite number of modes 
was proposed exploiting techniques from continuous variables. The covariance matrix 
formalism was employed to study the Unruh effect and extraction of entanglement from 
quantum fields without perturbation theory. The techniques introduced in [12] are 
restricted to simple trajectories in which the time evolution is trivial. To show in 
detail how the formalism introduced was applied, the authors presented over-simplified 
examples using detectors coupled to a single mode of the field which is formally only 
applicable when the field can be decomposed into a discrete set of modes with large 
frequency separation. This situation occurs, for example, when the detectors are inside 
a cavity. The detector model introduced in this work generalizes the model presented 
in [12] to include situations in which the time evolution is non-trivial. 

We introduce the mathematical techniques required to solve the time evolution of a 
detector which couples to an arbitrary time-dependent frequency distribution of modes. 
The detector is assumed to be again a harmonic oscillator. The interaction of the 
detector with the field is purely quadratic in the operators and, therefore, we can employ 
the formalism of continuous variables taking advantage of the powerful mathematical 
techniques that have been developed in the past decade [2]. Our techniques allow us 
to obtain the explicit time dependent expectation value of relevant observables, such 
as mean excitation number of particles. As a concrete example, we employ our model 
to analyze the response of a detector, which moves along a arbitrary trajectory and is 
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coupled to a time-dependent frequency distribution of field modes. Recently it was 
shown that a spatially dependent coupling strength can be engineered to couple a 
detector to a gaussian distribution of frequency modes [13]. Here we analyze the case 
where the coupling strength varies in space and time in such way that the detectors 
effectively couples to a time evolving frequency distribution of plane waves that can 
be described by a single mode. This technique simphfies the Hamitonian and an exact 
expression for the number operators can be found. We also discuss the extent of the 
impact of the techniques developed in this paper: in particular, we stress that they 
can be successfully applied for a finite number detectors following arbitrary trajectories. 
The formalism is also applicable when the detectors are confined within cavities. In 
this last case, the complexity of our techniques further simplifies due to the discrete 
structure of the energy spectrum. Finally, we note that the model can be generalized 
to the case where the detector is a quantum field itself. 

2. Interacting systems for relativistic quantum information processing 

Work in the field of relativistic quantum information has typically considered the 
standard point-like Unruh-DeWitt interaction between a detector and a field [10, 14, 15]. 
In this scenario the detector interacts with equal strength with all the modes of the field. 
In general, the exact equations that govern the dynamics of the system arc difficult to 
solve. Therefore, perturbation theory techniques are employed to evolve the system. 
Here we present a novel method to treat the exact dynamics of the interaction using 
methods from symplectic geometry. 

The general interaction Hamiltonian Hi{t) between a quantum mechanical system 
(detector) interacting with a bosonic quantum field x) in 4-dimensional spacetime 
is 



where (r, x) are a suitable choice of coordinates for this spacetime, m(r) is the monopole 
moment of the detector and g denotes the determinant of the metric tensor [16]. The 
function J-'(r, x) is the effective interaction strength between the detector and the field. 
When written in momentum space, it describes how the internal degrees of freedom 
of the detector couple to a time dependant distribution of the field modes. A specific 
example will be given in section 5. Such details can be determined by a particular 
physical model of interest. A more detailed and complete treatment of the interaction 
Hamiltonian (1) can be found in [13, 17, 18]. 

We consider a bosonic quantum field $ which can be expanded in terms of a 
particular set of solutions to the field equations 0k (t, x) as 



k 

where the variable k is a set of discrete parameters and Dy_ are bosonic operators that 
satisfy the canonical commutation relations [Dk^-Oj^./] = 5kk'- We refer to the solutions 
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0k as field modes. We emphasize that the modes 0k need not be standard solutions to 
the field equations (i.e. plane waves in the case of a scalar field in Minkowski spacetime) 
but can also be wave-packets formed by linear superpositions of plane waves. The modes 
are all normalized according to the appropriate inner product [19]. 
We can engineer the function F(t, x) such that 

j dVF{T, x)$(r, x) = /i(T)L>k. + h.c, (3) 

where one mode, labelled via k*, has been selected out of the set {0k}, which in turn 
implies 

Hi{r)=m{r)[h{r)D^,+li.c.]. (4) 

Therefore, the coupling strength has been specially designed to make the detector couple 
to a single mode. In the case of a free massless 1 + 1-dimensional scalar field, the mode 
the detector couples to corresponds to a time dependent frequency distributions of 
plane waves. More details about this decomposition and how to interpret the modes as 
frequency distributions are available in [4] . Note that for the case of a field contained 
within a cavity, where the modes are explicitly discrete, our methods apply in a simple 
manner. 

In the following, we explain how to derive the time evolution of an arbitrary 
number of detectors interacting with an arbitrary number of fields when the interaction 
Hamiltonian is of a purely quadratic form given by equation 4. 



3. Time evolution of N interacting bosonic systems 

The generalization of the quadratic Hamiltonian given by equation (4) to N interacting 
bosons is 

N{2N+l) 

H{t)= Y. ^At)G„ (5) 

where the functions Xj are real and the operators Gj are Hermitian and quadratic in 
creation and annihilation operators {{dj,dj)}. For example, G — d\dl + did2- The 
operators Gi form a closed Lie algebra with Lie bracket [Gj, Gj] — €ijkGk- The complex 
numbers eijk are called structure constants and form a completely antisymmetric matrix. 

We wish to find the time evolution of our interacting system. In the general case, 
the Hamiltonian H(t) does not commute with itself at different times [H{t),H{t')] ^ 0. 
Therefore, the time evolution is induced by the unitary operator 

C/(i) = ^e-'/o'^*W) (6) 

where ^ stands for time ordering. We can employ techniques from Lie algebra and 
symplectic geometry [20,21] to explicitly find a solution to equation (6). The unitary 
evolution of the Hamiltonian can be written as [22], 

U{t)^llU,{t)^lle-^^^('^^^ (7) 
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where the functions Fj[t) associated with generators Gj are real and depend on time. 
By equating (6) and (7) and differentiating with respect to time we obtain 

H{t)U{t) = Fid n Ujit) + e-^^^(*)«^F2(i)G2 \{ Uj{t) + ... (8) 

j i>2 

Multiplying on the right by U~^{t) gives us a sum of similarity transformations 

H(t) = F,(t)Gi + F2{t)UiG2Ur' + F^.(t)UiU2GsU^^Ur^ + ... (9) 

We have obtained a set of A^(2A^+1) coupled, non-linear, first order ordinary differential 
equations of the form 



j j 



^^At)F,{t) + V A,F;(t) + 7^{t) = 0. (10) 



The aik{t) and /3ik{t) can be functions of Fj{t) and Xj{t). The Hamiltonian and the 
initial conditions Fj{0) — completely determine the unitary time evolution operator 
(6). 

The equations can be re-written in a formalism which simplifies calculations by 
defining a vector that collects bosonic operators 

X:= (o?i,4,...,div,<)'^ (11) 

In this formalism the successive applications of the Baker-Campbell-Hausdorff formula 
which are required in the similarity transformations will be replaced by simple matrix 
multiplications reducing the problem from a tedious Hilbert space computation to simple 
hnear algebra. We write 

Uj{t) Gk U-\t) = • Sj(i)^ • Gj ■ Sj{t) ■ X (12) 

where 

Uj{t)^Uj\t) := Sj(t) -X 

and Gj is the matrix representation of Gj defined via Gj := X''" • Gj • X. The dynamical 
transformation of the vector of operators X generated by the interaction Hamiltonian 
Gj is given by the symplectic matrix [23] 

where Fj{t) are real functions associated with the generator Gj and Vtij := [Xj,Xj] is 
the symplectic form. The symplectic matrix S satisfies 8^08 = 0. In this formalism 
equation (9) takes the form 

X^ • H • X = Fi{t)X^ • Gi • X + F2(t)X^ • Si{ty • G2 • 8i(i) • X + . . . 

and therefore, the matrix representation of the Hamiltonian H is given by 

H = Fi(i)Gi + F2(i)Si(i)^ • G2 • Si(i) 

+ Fs{t)Siity ■ S2ity ■ G3 • 82(t) ■Si{t) + ... (14) 
The matrix products of the form 

Fj{t)Sk{ty -GySkit) (15) 
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in equation (14) must be explicitly computed such that it is possible to re- write equation 
(14) in terms of the generators Gj. By equating the coefficients of equation (14) to 
the coefficients A-,(r) in equation (5) we obtain a set of coupled N{2N + 1) ordinary 
differential equations. Solving for the functions Fj{t), we obtain the explicit expression 
for the time evolution of the system as described by equation (7). The final expression 
is 

S(t) = Sj{t) = Yl e-^^.(*)«G,, (16) 
j j 

which corresponds to the time evolution of the whole system. A case of great interest is 
that of Gaussian states which are common in quantum optics and relativistic quantum 
theory [2]. In this case the state of the system is encoded by the first moments (X^) 
and a covariance matrix T{t) defined by 

r,,- = (X,X,- + X,X,) - 2(X,)(X,). (17) 

The time evolution of the state is given by 

T{t) = sHt)mm- (18) 

4. Application: Time evolution of a detector coupled to a field 

We now apply our formalism to describe a situation of great interest to the field of 
relativistic quantum information: a single detector following a general trajectory and 
interacting with quantum field via a general time and space dependent coupling strength. 
For simplicity, consider an uncharged massless scalar field ^{t,x) in 1 + 1 dimensional 
Minkowski spacetime with coordinates {t,x) and metric g^i, — (—1, 1). The field obeys 
the standard Klein-Gordon equation = where □ := —df + d^ is the d'Alambertian 
in fiat Minkowski coordinates. The plane- wave solutions to the field equation are 

0,(t,x) = -^e--(*--) (19) 

which are (Dirac delta) normalized {(f)^^, 0^/) = 6{u — u') through the standard conserved 
inner product (■, ■), see [19]. The field can be expanded in terms of these solutions as 

r+oo 

/ da;M,(i,x) + at,C(i,x)], (20) 
Jo 

The annihilation operators a^^ define the well known Minkowski vacuum |0)^ through 

ajO);^ = 0, > 0. (21) 

The degrees of freedom of the detector which we assume to be a harmonic oscillator are 
described by the bosonic operators d, d^ that satisfy the usual commutation relations 
[d^d^ = 1. The vacuum |0)^ of the detector is defined by d\Q)^ = 0. Therefore, the 
vacuum |0) of the non-interacting theory takes the form |0) := lO)^^ (8) \0)m- 

In the interaction picture, we assume that the detector couples to the field via the 
interaction Hamiltonian 

/r+oo 
dx^T{t, x) dw [aMt: ^) + ^)] > (22) 
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where x) is a time and space dependent coupling strength and g denotes the 
determinant of the metric tensor. In momentum space we can see that the detector 
couples to a time dependent frequency distribution of planes waves 

Hi{t) = m{t) J dujg{t,u)e''^* (23) 

where 

g{T,u) ^ j dx^f^F{T,x)e-''''' (24) 

is the frequency window of the detector which in general evolves in time. We now 
parameterise the interaction via a suitable set of coordinates, (r, ^), that describe a 
frame comoving with the detector [17, 18]. 

In this comoving frame, the monopole moment of the detector is 

m(T) = s{t) [e-'^^d + e'^^d^ (26) 

s(t) is real function which allows us to turn on and off the interaction and A is 
frequency of the harmonic oscillator. In momentum space the detector couples to a 
time-dependent frequency distribution of Minkowksi plane-wave field modes. In [13] the 
spatial dependence of the coupling strength was specially designed to couple the detector 
to peaked distributions of Minkowsky or Rindler modes. Here we consider a coupling 
strength that can be designed to couple the detector to a time- varying wave-packet. It 
is therefore more convenient to decompose the field not in the plane-wave basis but in 
a special decomposition 

$ = + Dhh* + (27) 

where ip is the mode the detector couples to, which corresponds to a time dependent 
frequency distribution of plane waves. The operators D, D'^ arc time independent and 
satisfy the canonical commutation relations [D^D'^] = 1. The field $' includes all the 
modes orthogonal to 'ip and we will assume them to be countable. Therefore 

h{r) = j d^J^{r,0^{r,0, j d^Hr^i)^' {t,^) = Vr. (28) 

Such a decomposition can always be formed from a complete orthonormal basis (an 
example of which can be found in [17]). In general, the operator D does not annihilate 
the Minkowski vacuum |0)jy^. This observation is a consequence of fundamental ideas 
that lie at the foundation of quantum field theory, where different and inequivalent 
definitions of particles can coexist. Such concepts are, for example, at the very core of 
the Unruh effect [10]. The operator D will annihilate the vacuum |0)^. Note that the 
vacuum state |0)y^ of this interacting system is different from the vacuum state |0) of the 
noninteracting theory, i.e. |0) ^ |0)j. 

In this decomposition, the interaction Hamiltonian takes a very simple form 

Hj{t) = m(r) ■ [h{T)D + h*{T)D^] (29) 



Time evolution techniques for detectors in relativistic quantum information 



8 



which describes the effective interaction between the internal degrees of freedom of 
a detector following a general trajectory and coupling to a single mode of the field 
described by D. The time evolution of the system can be solved in this case by employing 
the techniques we introduced in the previous section. However, this formalism is directly 
applicable to describe the interaction of N detectors with the field. In that case, our 
techniques yield differential equations which can be solved numerically. We choose 
here to demonstrate our techniques with the single detector case since it is possible to 
compute a simple expression to the expectation value of the number of particles in the 
detector. 

Let the detector-field system be in the ground state |0)^ at r = 0. We design a 
coupling such that we obtain an interaction of the form (29). In this case, the covariance 
matrix only changes for the detector and our preferred mode. The subsystem described 
by D is always separable from the rest of the non-interacting modes. The covariance 
matrix of the vacuum state |0)^ is represented by the 4x4 identity matrix, i.e. r(0) = 1. 
The final state T{t) therefore takes the simple form of T{t) = S^S. The final state 
provides the information we need to compute the time dependent expectation value 



Na{T) := {Sd){T). 

Prom the definition of the covariance matrix T{t), one finds that Nci{t) is related 
to T{t) by 



In this paper we choose to work with states that have first moments zero, i.e {Xj) — 0. 
In this case, since our interaction is quadratic, the first moments will remain zero [24]. 
Therefore we are left with equation 



Our expressions hold for detectors moving along an arbitrary trajectory and coupled 
to an arbitrary wave-packet. Given a scenario of interest, one can solve the differential 
equations, obtain the functions Fj(t) and, by using the decomposition in equation (16), 
one can obtain the time evolution of the system. We can then find the expression for 
the average number of excitations in the detector at time r, which reads 



For this choice of initial state we find that the functions -Fj(r) are associated with the 



generators d = d^D^ + dD,G2 = -i{d^D^ - dD),G^ = d^^ + d^,G4 = -i{d^^ - d^), 



respectively. The appearance of these functions can be simply related to the physical 
significance of the operators Gj. In fact, the generators Gi and G2 are related to 
the well known two-mode squeezing operators. Such operations generate entanglement 
and are known to break particle number conservation. The two generators G3 and G4 
arc related to the single-mode squeezing operators for the mode d. The generators 
G1...G4, together with the generators G5 = D^'^ + D"^ and Gq = —i{D^'^ — D"^) 
which represent single mode squeezing for the mode form the set of so-called active 
transformations of a gaussian state and do not conserve total particle number. The 




(30) 




(31) 



Nd{T) 



cosh(2Fi(T))cosh(2F2(T))cosh(2F3(T))cosh(2F4(T)) - 1 

2 



(32) 
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remaining operators, whose functions are absent in equation (32), form the passive 
transformations for gaussian states. These operators are also known as generahsed 
beam sphtter [23] ; they conserve the total particle number of a state and hence do not 
contribute to equation (32). 

It is also of great physical interest to study the response of one (or more) detector 
when the initial state is a different vacua |0), for example the Minkowski vacuum |0)m- 
It is well known [16] that the relation between different vacua, for example |0) and \Q)d) 
is 

|0) = Are-^^^^-^*^-°'^-^l|0)^, (33) 

where A?" is a normalization constant and the symmetric matrix V is related to the 
Bogoliubov transformations between two different set of modes {(l>uj}-i{'4^j} used in the 
expansion of the field. The modes {0u,} are related to the vacuum |0) while the 
modes {■0^} are related to the vacuum |0)^. In general, the matrix V takes the form 
V := B*A~^ [16], where the matrices A and B collect the Bogoliubov coefficients 
^joj) -Sja; which for uncharged scalar fields are defined as 

^,-u, = (V',-,0j (34) 
S,.^ = -(V',-,C) (35) 
In the covariance matrix formalism, we collect the detector operators d, S and mode 
operators Dj, D\ in the vector X := (d, Di, dI, D2, D\, . . .f. The initial state r(0) 
will not be the identity anymore, r(0) — A 1. The final state T{t) will take the form 

r(T) = St(T)AS(T) (36) 

We assume that within the field expansion using the basis {ipj}, the detector interacts 
only with the mode ipp. We can again apply our techniques to obtain the number 
expectation value as 

Na{t) = ^^^12 + ^!(2+,)^i(2+p) (VVt)^^ 

~ ^i{2+p)^H3+p) V*p - 'S'J(3_^p)5'i(2+p) Vpp. (37) 

5. Concrete Example: inertial detector interacting with a time-dependent 
wavepacket 

To further specify our example we consider the detector stationary and interacting with 
a localized time dependent frequency distribution of plane waves. The free scalar field 
is decomposed into wave packets of the form [17] 

u^:=fdkf^(kW (38) 

where the distributions fmi{k) are defined as 

J n^-s ._ ; e-i/2e-2-'^A e(m - 1/2) <k<e{m + 1/2) 
' otherwise 
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with e > and {m, /} running over all integers. The mode operators associated with 
these modes are defined as 

ami := j dkf^i{k)ak (40) 

The distributions fmi{k) satisfy the completeness and orthogonality relations 

E fmi(k)fUk') = S{k - k') j dkfmi{k)f*r^,r{k) = 5mm'5iv (41) 



The wave packets are normalised as {umhUm'i') = ^mm'^w and operators satisfy the 
commutation relations [o-ni;- = 5,, 
terms of these wave packets as 

^ = [^miami + h.c 

ml 

Following [18], we consider an inertial detector and so we can parametrise our interaction 
via t — T and x = ^. We now construct our detectors spatial profile to be 



mm'^w- The scalar field can then be expanded in 

(42) 



g{T,i) h{T) j dkf*^Muk{r,0 



(43) 



where /i(t) is now an arbitrary time dependent function which dictates when to switch 
on and off the detector. Physically, this corresponds to a detector interaction strength 
that is changing in time to match our preferred mode, labelled by M, L. We point 
out that any other wave packet decomposition could be chosen as long as it satisfies 
completeness and orthogonality relations of the form (41). The form of the spatial 
profile to pick out these modes is therefore general. Inserting the profile (43) into our 
interaction Hamiltonian (25), we obtain 

Hj{t) = {de-'^^ + Se+'^^) (h{r)aML + h*{r)alr) (44) 

We choose the switching on function to be /i(t) = r^e""^^/^, where T modulates the 
interaction time. The interaction Hamiltonian is then 

Ht{t) = T^e-^'l^ {de-'^^ + rfte+^^-) (a^M + a[^) (45) 

which written in symplectic form (see equation (5)) is 

Hi{t) = r'e-^'/^ [cos (tA) d + sin (tA) G2 

+ cos (tA) Gr + sin (tA) Gg] (46) 

where = d^D + dD^ and Gs — —i{dD^ — d^D). The matrix representation of Hj is 











irA 
itA 



g-irA g«T"A 













(47) 



Equating (47) and (14), or equivalently (46) and (5), gives us the ordinary differential 
equations we need to find the functions Fj for this specific example. Here we solve the 
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tau 

Figure 1. Mean number of particles, ndir), as a function of time t. Here we used 
T = 80 and A = 2tt. 

equations for Fj{t) numerically and we plot the average number of detector excitations 
nrf(r) as a function of time in figure (5). 

We find that the number expectation value of the detector grows and oscillates as 
a function of time while the detector and field are coupled. This can be expected since 
the time dependence of the Hamiltonian is multiplied by complex exponentials that 
will induce phase rotations in the state and hence oscillations in the number operator. 
Finally, the number expectation value reaches a constant value after the interaction 
is turned off. This is expected as well since, once the detector evolves freely, the free 
Hamiltonian does not account for emissions of particles from the detector. 

6. Discussion 

It is of great interest to solve the time evolution of interacting bosonic quantum systems 
since they are relevant to quantum optics, quantum field theory and relativistic quantum 
information, among many research fields. In most cases, it is necessary to employ 
perturbative techniques which assume a weak coupling between the bosonic systems. In 
relativistic quantum information, perturbative calculations used to study tasks such as 
teleportation [11] and extraction of vacuum entanglement [11] become very complicated 
already for two or three detectors interacting with a quantum field. In cases where the 
computations become involved, physically motivated or ad hoc approximations can aid, 
however, in most cases, powerful numerical methods must be invoked and employed to 
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study the time evolution of quantities of interest. 

We have provided mathematical methods to derive the differential equations that 
govern the time evolution of N interacting bosonic modes coupled by a purely quadratic 
interaction. The techniques we introduce allow for the study of such systems beyond 
perturbative regimes. The number of coupled differential equations to solve is A^(2iV+l), 
therefore making the problem only polynomially hard. Symmetries, separable subsets of 
interacting systems, among other situations can further reduce the number of differential 
equations. 

The Hamiltonians our method is apphcable to include a large class of interactions. 
In this paper, as a simple example, we have applied our mathematical tools to analyze 
the time evolution of a single harmonic oscillator detector interacting with a quantum 
field. However, our techniques are readily applied to N detectors following any trajectory 
while interacting with a finite number of wave-packets through an arbitrary interaction 
strength x). Our techniques simplify greatly when the detectors are confined within 
a cavity where the field spectrum becomes discrete. The cavity scenario allows one 
to couple a detector to a single mode of the field in a time independent way as, in 
principle, no discrete mode decomposition needs to be enforced. Therefore, the single 
mode interaction Hamiltonian (4) can arise in a straight forward fashion. Inside a 
cavity, the oversimplified examples introduced in [12] where two harmonic oscillators 
couple to a single mode of the field arc well known to hold trivially. 

We have further specified our example to analyze the case of an incrtial detector 
interacting with a time-dependent wave-packet. We have showed how to engineer a 
coupling strength such that the interaction Hamiltonian can be descried by an effective 
single field mode. However, the field mode is not a plane-wave but a time dependent 
frequency distribution of plane waves. In this case we have solved the differential 
equations numerically and showed the number of detector excitations oscillates in time 
while the detector is on. 

Work in progress includes using these detectors to extract field entanglement and 
perform quantum information tasks. 

Note 

Near the completion of this work but before posting our results, we became aware of 
another group working independently along similar lines [25]. We agreed to post our 
results simultaneously. 

7. Acknowledgments 

We would like the thank Jorma Louko, Achim Kempf, Sara Tavares, Nicolai Priis and 
Jandu Dradouma for invaluable discussions and comments. I. F. acknowledges support 
fi-om EPSRC (CAF Grant No. EP/G00496X/2). 



Time evolution techniques for detectors in relativistic quantum information 



13 



References 

[1] Nielsen M A and Chuang I L 2000 Quantum Computation and Quantum Information (Cambridge, 
N.Y.) 

Alsing P M and Fuentes I 2012 Classical and Quantum Cravity 29 224001 
Bruschi D E, Fuentes I and Louko J 2012 Phys. Rev. D 85 061701 

Dragan A, Doukas J, Martin-Martinez E and Bruschi D E 2012 Localised projective measurement of 
a relativistic quantum field in non-incrtial frames {Preprint arXiv: 1203.0655vl [quant -ph]) 
Downes T, Ralph T and Walk N 2012 Quantum communication with an accelerated partner 

{Preprint arXiv:1203.2716vl [quant-ph]) 
Priis N, Bruschi D E, Louko J and Fuentes I 2012 Phys. Rev. D 85 081701 
Friis N, Huber M, Fuentes I and Bruschi D E 2012 Phys. Rev. D 86 105003 

Bruschi D E, Dragan A, Lee A R, Fuentes I and Louko J 2012 Motion-generated quantum gates 

and entanglement resonance {Preprint arXiv: 1201 .0663v2 [quant-ph]) 
Bruschi D E, Louko J, Faccio D and Fuentes I 2012 Mode-mixing quantum gates and entanglement 
without particle creation in periodically accelerated cavities {Preprint arXiv: 1210.6772 
[quant-ph] ) 
Unruh W G Phys. Rev. D 14 870 

Lin S Y, Kazutomu S, Chou C H and Hu B L 2012 Quantum teleportation between moving 

detectors in a quantum field {Preprint 1204.1525 [quant-ph]) 
Dragan A and Fuentes I 2012 Probing the spacetime structure of vacuum entanglement {Preprint 

arXiv: 1105. 1192 [quant-ph]) 
Lec A R and Fuentes I 2012 Spatially extended unruh-dewitt detectors for relativistic quantum 

information {Preprint arXiv: 1211 . 5261 [quant-ph]) 
DeWitt B S 1979 Quantum gravity: The new synthesis General Relativity: An Einstein Centenary 

Survey ed Hawking S W and Israel W (CUP Archive) 
Grove P G and Ottewill A C 1983 Journal of Physics A: Mathematical and General 16 3905 
Fabbri A and Navarro-Salas J 2005 Modeling Black Hole Evaporation (Imperial College Press) 
Takagi S 1986 Progress of Theoretical Physics Supplement 88 1-142 
Schlicht S 2004 Class. Quantum. Grav. 21 

Crispino L C B, Higuchi A and Matsas G E A 2008 Rev. Mod. Phys. 80 
Berndt R 2001 An Introduction to Symplectic Geometry (American Mathematical Soc.) 
Hall B 2003 Lie Groups, Lie Algebras, and Representations: An Elementary Introduction 
(Springer) 

Puri R R 2001 Mathematical Methods of Quantum Optics (Springer) 

Luis A and Sanchez-Soto L L 1995 Quantum and Semiclassical Optics: Journal of the European 

Optical Society Part B 7 
Weedbrook C, Pirandola S, Garcia-Patron R, Cerf N J, Ralph T C, Shapiro J H and Lloyd S 2012 
Rev. Mod. Phys. 84 

Brown E G, Martin-Martinez E, Menicucci N C and Mann R B 2012 Detectors for 
probing relativistic quantum physics beyond perturbation theory {Preprint arXiv: 1212. 1973 
[quant-ph] ) 



